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We derive the explicit expression for the distribution of 
resonance widths in a chaotic quantum system coupled to con- 
tinua via M equivalent open channels. It describes a crossover 
from the x 2 distribution (regime of isolated resonances) to a 
broad power-like distribution typical for the regime of over- 
lapping resonances. The first moment is found to reproduce 
exactly the Moldauer-Simonius relation between the mean 
resonance width and the transmission coefficient. This fact 
may serve as another manifestation of equivalence between 
the spectral and the ensemble averaging. 

PACS numbers: 05.45. +b, 24.30 v 



Chaotic scattering has been a subject of rather inten- 
sive research activity during the last decade, both the- 
oretically and experimentally, see reviews 0]- H. This 
phenomenon is encountered in a variety of physical sys- 
tems ranging from atoms, atomic nuclei and molecules to 
mesoscopic devices [Q and microwave cavities |^|| . 

One of the basic concepts in the domain of chaotic 
quantum scattering is the notion of resonances, repre- 
senting long-lived intermediate states to which bound 
states of a "closed" system are converted due to cou- 
pling to continua. On a formal level resonances show 
up as poles of the scattering matrix S a b(E) occuring at 
complex energies Ek = £k — §1% where £t and Yk are 
called position and width of the resonance, correspond- 
ingly. The general problem of determining the domain 
of concentration and distribution of S-matrix poles is of 
fundamental interest in the scattering theory and is a 
field of intensive research activity Q| . Powerful numer- 
ical methods are available ( e.g the method of complex 
scaling) allowing one to extract resonance parameters for 
models in atomic and molecular physics B. 

Whereas the issue of energy level statistics in closed 
chaotic systems was addressed in an enormous amount 
of papers (see, e.g. || and references therein), statistical 
characteristics of resonances are much less studied and 
attracted significant attention only recently. In the case 
of weak effective coupling to continua individual reso- 
nances do not overlap: (r) <C A, with A standing for the 
mean level spacing of the " closed" system and (r) stand- 
ing for the mean resonance width. Under these conditions 
one can use a simple first order perturbation theory to 
calculate resonance widths in terms of eigenfunctions of 
the closed system, see e.g [fljlQl. Quite generally, one 
finds in such a procedure that the scaled widths y s = j^j 



are distributed according to the so-called ^-distribution: 

(1) 



I>/2) 

where T(z) stands for the Gamma function and the pa- 
rameter v = M (y — 2M) for systems with preserved 
(broken) time reversal symmetry (TRS), M being the 
number of open scattering channels. The basic fact un- 
derlying this simple form of eq. (|l|) is the Gaussian statis- 
tics of components of wavefunctions. The case v = 1 is 
known as Porter-Thomas distribution [jllj . It was shown 
to be in agreement with experimental data in neutron- 
nuclei resonances, the fluorescence excitation spectrum 
of the NO2 molecule, the diamagnetic Rydberg spectrum 
of the lithium atom (see references in M, Jl^ |, |11|) and, 
more recently in microwave cavities [jfj. 

When the coupling to continua increases some reso- 
nances start to overlap and the simple perturbative re- 
sult eq.(P loses its validity. This regime is most easily 
reached when the number of open channels M is of the 
same order as the number N of bound states in the closed 
system, the situation encountered in the so-called " three 
disk" scattering problem |Q. The range of parameters 
M oc N was studied in much detail in It corre- 

sponds to very strong overlap of resonances (r) >> A 
which is the opposite limiting case as compared to the 
domain of validity of the distribution eq.(Q). 

In the general case of M open channels, M very un- 
stable states (broad resonances) are formed, whereas the 
rest, N — M states, have much smaller widths ]l4] , p^[ . 
Changing the effective coupling to continua one arrives 
at the situation when a substantial portion of these long- 
lived resonances show a considerable overlap, even for 
the case of few open channels M <C N. Experimentally, 
one quite frequently encounter the case of only M ~ 1 
open channels and (r) ~ A PJip jlT^ , ^6[ | . In this situa- 
tion neither application of eq.W) nor exploitation of the 
results from |l3[| can be reasonably justified. The main 
goal of the present publication is to fill in this gap pro- 
viding the explicit distribution covering the whole range 
of the parameter (T)/A for the few channel scattering 
problem. 

At present, there are two complementary theoretical 
tools employed to calculate statistical characteristics of 
"open" quantum systems whose closed classical coun- 
terparts demonstrate chaotic behaviour. These are the 
semiclassical approach [|l]j^] and the stochastic approach 
]l7t , the relation between both methods being compre- 
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hensively discussed in In particular, the semiclassi- 
cal approach operates with genuine microscopic Hamilto- 
nians and allows for treating particular systems with full 
account of their specific features. However, the applica- 
bility of the method is essentially restricted to the case 
of many open channels M ^> 1. In contrast, the stochas- 
tic approach could, in principle, deal with as many open 
channels as one likes. As a price, it is able to address 
only universal, i.e. system independent characteristics 
common to most of quantum chaotic systems and stem- 
ming from the very basic fact of chaotic internal motion. 
Such kind of generality is achieved by incorporating the 
random matrix description of chaotic systems as an input 
information. The latter principle is a commonly accepted 
one in the domain of Quantum Chaos ||. 

The starting point of the stochastic approach is the 
following multichannel Breit-Wigner representation for 
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S ab (E) = 6 ab - 2inJ2w* a [E - TUftfWjt 

{n e f) ij = H i --mY Ja w ia w; a 



(2) 



Here the N x N matrix Hij is to model the Hamiltonian 
of the closed chaotic system and thus is chosen to be a 
member of the corresponding Gaussian ensemble of ran- 
dom matrices. To be specific, we consider in the present 
publication only systems with broken TRS by taking 
to be a random Hermitian matrix from the Gaussian Uni- 
tary Ensemble ||. The amplitudes W a i, a= 1, 2, M 
are matrix elements coupling the internal motion to one 
out of M open channels. Without much loss of generality 
these amplitudes can be chosen in a way ensuring that 
the average S— matrix is diagonal in the channel basis: 
Sab = 8 a bSaa- Provided the energy E is real, one finds 
the following expression: p7[: 



— _ 1 - lc g(E) 

baa -i + lc9 (Ey 
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where g(E) = iE/2+ (1 -E 2 /^) 1 / 2 and we assumed that 
all M channels are statistically equivalent for the sake of 
simplicity. The strength of coupling to continua is con- 
venient to characterize via the transmission coefficients 
T a = 1 — \S aa \ 2 that are given for the present case by the 
following expression: 



1 
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2Re g(E) 



(7c + 7 c l ) 
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The quantity T a measures the part of the flux in channel 
a that spends substantial part of the time in the inter- 
action region 0,^8). We see that both limits 7 C — > 
and 7 C — > oo equally correspond to the weak effective 
coupling regime T a < 1 whereas the strongest coupling 
(at fixed energy E ) corresponds to the value -f c = 1. 
The maximal possible coupling corresponding to the up- 
per bound T a — 1 is achieved in the present model for an 
energy interval in the vicinity of the center E = 0. 



In order to get access to the distribution of resonance 
widths we use the fact that resonances are actually eigen- 
values of the non- Hermitian Hamiltonian matrix TL e f de- 
fined in eq.(|| ). Two-dimensional density of these eigen- 
values p(X, Y) in the complex plane E = X + iY can be 
found if one knows the "potential" fl3|| : 



lnDet[(E - H ef )(E - W e/ )+ + k 2 } 



2ttN 



in view of the relation: p(X,Y) — lim K ^o d 2 Y, n), 
where d 2 stands for the two-dimensional Laplacian. 
Technically, it turns out to be much easier to restore the 
potential from its derivative: 



~ 1 d lim — lnZ( Kb k) 
dn 2 2nN dn K b ^K Ok 

7Ik k \ _ Det[(E-H sf )(E-H <!f y+ K 2 ] 
A » 6 ' K > ~ Det[(E-H cf )(E-H cf Y>+ K i] 



(5) 



Then one can find the following representation for the 
generating function Z(k , k) in terms of the Gaussian in- 
tegral over both commuting and anticommuting (Grass- 
mann) variables: 

(-l) N Z( Kb , k) = /[d*] exp{-C Q (^) - AW} 

£i(*) = -Y(^& y) - *t(f <g) ct )* + (k- K b )(&K$) 

(6) 

where 

*t = JslJlxlxl)\ m = Y[ p=ia dS P dSld Xl d x i 
with S p and \v being N— component vectors of complex 
commuting and Grassmannian variables, respectively. 
The 4x4 matrices A, L, &o and K are (block)diagonal 
of the following structure: 



A = diag(l, -1,1,-1); 
K = diag(0, 0, 1, -1); 



L = diag(l, -1,1,1) 
(fo = diag(iS :c ,S a; ) 



and Yi x = 



The entries of the N x N matrix f 
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1 / 

are given by = -it J2 a w ^ w * a . 

It is easy to notice that if one suppresses the second 
term £i( v I') in the exponent of eq.(^|), the corresponding 
Gaussian (super)integral just coincides with the gener- 
ating function emerging in calculations of the two-level 
spectral correlation function for the pure GUE |l9| . This 
fact allows one to apply the well- developed technique 
pioneered by Efetov | 19| in order to perform the av- 
eraging over the GUE of matrices H . After a set of 
standard manipulations |l7],|l9| one arrives at the fol- 
lowing expression for the density p(y) of level widths 
y = ^ = — 2tt~p(X) NY of those resonances, whose posi- 
tions are within the narrow window around the point X 
of the spectrum pfj : 



p(y) 



1 d 2 
2^V 
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(7) 
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where 



' dz 



<i>{y,z) = 

dn(Q)Str(KQ)Sder M / 2 

x exp{-f Str(QA) - l -§Str(Qa)} 



A 2vr t 



1 M 

_ = -_ln(l-T B (10) 

1 Z7T 



(8) 



where ct = diag(E;j, E x ) and the (graded) matrices Q sat- 
isfying Q 2 — — 1 are taken from the graded coset space 
[/(l,l/2)/[/(l/l)(g) [7(1/ , 1 ) , whose explicit par ametriza- 
tion can be found in JL9| . Here we introduce the notation 
){lc + TcT^and used the symbols Str,Sdet 



' ~ 2Trp(X) ' 

for the graded trace and the graded determinant, corre- 
spondingly. 

Now it is evident, that the form of the resonance widths 
distribution does not depend on S a b itself. Rather, it 
is uniquely determined by the value of the transmission 
coefficient equal to T a = 2/(1 + r) for the present case. 
This is very satisfactory feature shared by all observables 
like crossections, etc. and believed to be an indicator of 
universality of the obtained results jl7|,[l8) ■ 

Still, the evaluation of the superintegral in eq.(^|) and 
subsequent restoration of the density p{y) is quite an 
elaborate task. For this reason we choose to skip all 
technical details in favour of presenting the outcoming 
result, which turns out to be surprisingly simple: 



p(v = £ ) = 



f T iV dtt 

yr sinh y 



M, 
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This expression provides the explicit form for the dis- 
tribution of level widths for a M-channel chaotic system 
with broken time-reversal invariance. It constitutes the 
main result of the present publication. 

Quick inspection of eq.(^|) shows that it is indeed re- 
duced to the x 2 distribution, eq.(|l]) when the effective 
coupling to continua is weak: t»1. Under this condi- 
tion resonances are typically too narrow to overlap with 
others: y < t^ 1 <C 1. However, as long as the effective 
coupling becomes stronger, the parameter r decreases 
towards unity. Under these conditions another domain 
of resonance widths becomes more and more important: 
7q~L < V < 7ZT> where the distribution eq.(^) shows the 
powerlaw decrease: p(y) ~ \My~ 2 . The most drastic 
difference from eq. (|l|) occurs for the maximal effective 
coupling t = 1 (i.e j c = 1 and E = 0). In this regime the 
powerlaw tail extends up to infinity making all positive 
moments ( starting from the first one) to be apparently 
divergent Q. 

It is interesting to note that the behaviour p(y) ss 
\My~ 2 nicely matches that obtained in |l3| for the 
asymptotic limit M oc JV > 1. Another point which 
is worth to be mentioned is that evaluating the first mo- 
ment of the distribution eq.(^|) exactly we arrive at the 
following expression for the mean resonance widths: 



This formula is well known in nuclear physics as 
Moldauer-Simonius relation plfl . It was derived for sys- 
tems with preserved TRS by averaging the S— matrix 
over the energy spectrum and using the unitarity con- 
dition. The fact that the ensemble averaging produces 
exactly the same expression is an explicit manifestation 
of the validity of the ergodicity hypothesis and can be 
considered as another support for the applicability of the 
stochastic approach. 

Generally, we see, that the logarithmic divergency of 
(r) at the critical coupling r = 1 is a direct consequence 
of the 1/y 2 decrease of the probability distribution eq.(|]). 
Thus, one may say, that such a powerlaw tail is typical 
for chaotic systems strongly coupled to continua and is 
just one of the clear manifestations of the strong overlap 
between individual resonances. 

To this end, it is interesting to mention that a little dif- 
ferent powerlaw distribution of resonance widths - that of 
the form p(y) cx y~ 3 / 2 - was observed in numerical stud- 
ies of the "diffusive " quantum chaotic system strongly 
coupled to continua |Q . The authors suggested a trans- 
parent qualitative explanation of this effect based on the 
very fact of classical diffusive dynamics. This argumen- 
tation, when appropriately modified, allows one to un- 
derstand the form of the powerlaw tail in our case as 
well. Indeed, it is well-known that the choice of GUE as 
a model for the closed chaotic system precludes the dif- 
fusive spread of wavepackets from being taken into con- 
sideration. Instead, the corresponding "chaotic" classical 
dynamics can be visualised as that of a particle hopping 
in a phase space of N states. The one-step transition 
probabilities between any pair of these states are equal. 
Among the totality of N states there are M states strongly 
coupled to continua: the particle is " absorbed" (escapes 
from the sample) when it hits any of these M states. 
Assuming M/N <C 1 the probability to survive up to a 
time r is simply p(r) oc exp [-|t]. Then the distribu- 
tion of "inverse escape time" r = r _1 can be estimated 

as p(T) = dp ^ r ^ which immediately gives p(y) oc M/y 2 
for the normalized width y = NT, provided y ^> M. 

The following comment is appropriate here. Dealing 
with realistic models of chaotic systems containing no 
random parameters one always performs statistics over 
an energy interval 5E containing many levels: SE ^> A, 
but being small enough for a systematic variation of the 
smoothed level density p(E) to be neglected: SE <C 
~p~l One may expect that universal features of such 
a statistics are adequately reproduced within the frame- 
work of the stochastic approach, but only on the level 
of " local- in-spectrum" characteristics calculated at fixed 
value of E. Indeed, any spectral averaging in the stochas- 
tic model performed on a scale comparable with the ra- 
dius of the semicircle unavoidingly mixes up data corre- 
sponding to very different values of the transmission coef- 
ficients, the procedure washing out any relevant physical 
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information. In particular, it seems quite meaningless to 
consider quantities like the " globally" averaged resonance 
width r g f = -k ^2 k Tfc which can be trivially found from 
the sum rule: T g i/A = j^Tr lmH e f = 2Mj c . This 
quantity can not be related to any particular transmis- 
sion coefficient. This fact, however, should not be mis- 
interpreted as impossibility to have universal statistics 
of S— matrix poles within the stochastic approach Jl8| . 
Rather, the quantity T g i can be found via the direct in- 
tegration of the universal expression (F)/A , eq. (|l0|) , over 
the energy E, upon substituting there the energy depen- 
dent value T a (E) from eq.(Q). The result of the integra- 
tion is finite for any j c , thus concealing a specific role 
of the critical coupling "f c = 1 when resonances with a 
divergent local mean width occur sufficiently close to the 
center of the spectrum |24|] 

The best candidates for checking the applicability of 
eq. (|]) to real physical systems are realistic models of bal- 
listic mesoscopic devices subject to an applied magnetic 
field that serves to break the TRS [Q . Very recently there 
were also reports on realizing TRS breaking in microwave 
cavities p5| . However, it is not sufficiently clear whether 
the specific ways to violate TRS used in the cited papers 
can be simulated by simple substitution of the GUE ma- 
trices for H in the framework of the stochastic approach. 

The best experimentally accessible resonance data @ 
as well as models used to study chaotic scattering corre- 
spond to preserved TRS. In order to treat this situation 
in the framework of the stochastic approach one has to 
use GOE rather than GUE matrices. All basic steps lead- 
ing to eq. (0) remain essentially the same, but the graded 
matrices Q acquire larger size making the analytical com- 
putations very involved. It is however quite clear that all 
the basic qualitative features of the distribution eq.(||) 
(in particular, the powerlaw behaviour p(y) oc My~ 2 for 
the overlapping resonance regime) should show up in that 
case as well. This fact is already indicated by the valid- 
ity of the Moldauer-Simonius relation for such systems. 
Moreover, this relation gives support to a conjecture, that 
the expression eq. (^) may well turn out to be a good ap- 
proximation for M— channel systems with preserved TRS 
provided one substitutes M/2 for M into the first form 
of the eq. (^) and interprete y as y — |J ■ 
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